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The catalysis of nuclear reactions by negatively charged relics leads to increased 
p^ . outputs of primordial ^Li and ^Be. In combination with observational constraints 

'^ \ on the primordial fractions of ^Li and ^Be, this imposes strong restrictions on the 

C^^ ' primordial abundance and the lifetime of charged relics. We analyze the constraints 

!r^ . from the catalysis of ^Be on supersymmetric models in which the gravitino is the 

Q \ lightest supersymmetric particle and a charged slepton — such as the lighter stau — the 

next-to-lightest supersymmetric particle (NLSP). Barring the special cases in which 
the primordial fraction of the slepton NLSP is significantly depleted, we find that 
the ^Be data require a slepton NLSP lifetime of less than 6 x 10'^ s. We also address 
C^ \ the issue of the catalytic destruction of ^Li and ^Be by late forming bound states of 

protons with negatively charged relics finding that it does not lead to any significant 
modification of the limit on the slepton lifetime. 
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1 Introduction 

Physics of the early Universe keeps proving to be an invaluable tool for probing particle 
physics models and, in particular, models of New Physics beyond the Standard Model. 
One of the most important viability checks of these models results from the epoch of Big 
Bang Nucleosynthesis (BBN), i.e., from the early Universe at temperatures of T < 1 MeV. 
The combination of Standard Model physics, general relativity, and the experimental de- 
termination of the baryon-to-photon ratio from the anisotropics of the cosmic microwave 
background [1] forms the framework for the standard BBN (SBBN) which makes firm pre- 
dictions for the primordial abundances of light elements such as deuterium, helium, and 
lithium. The comparison of SBBN predictions with observationally determined primordial 
fractions of these elements provides an important consistency check of standard cosmology 
and serves as a remarkably powerful discriminator among models of New Physics [2]. The 
three most generic ways in which New Physics can affect the outcome of BBN are the 
change in the timing of the reactions caused, e.g., by new significant contributions to the 
Hubble expansion rate [2], the non-thermal processes from late annihilation and decays of 
heavy particles [3, 4, 5, 6, 7, 8, 9, 10, 11], and the thermal catalysis of nuclear reactions 
caused by electromagnetically or strongly interacting relics [12]. 

Catalyzed Big Bang Nucleosynthesis (CBBN) was already discussed almost twenty years 
ago in Ref. [13, 14] . However, only in the last two years after the appearance of Ref. [12] , a lot 
of work has been done [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33] 
in order to refine various aspects of catalysis and to understand the implications of CBBN 
in the framework of specific models. The significant interest in CBBN is also fuelled by 
its direct connection with collider physics. Indeed, the relics causing catalysis could be 
directly produced at the Large Hadron Collider (LHC) because of their electromagnetic 
and/or strong interactions. One of the most interesting and perhaps one of the most natural 
frameworks is in this context the supersymmetric (SUSY) extension of the Standard Model. 
There a spectrum with the gravitino G as the lightest supersymmetric particle (LSP) and 
a charged scalar lepton li — such as the lighter stau ri — as the next-to-lightest (NLSP) is 
a commonplace occurence even if one adopts restrictive assumptions on the soft SUSY 
breaking sector [34, 35, 36, 17, 19, 26, 27, 28]. While the gravitino LSP is a promising 
candidate for dark matter [37, 38, 39, 40, 41, 42, 43, 44], the charged slepton NLSP can be 
long-lived and thus lead to CBBN [12, 17, 18, 19, 20, 21, 22, 23, 25, 26, 27, 28, 32]. If such a 
scenario is realized in nature, each Standard Model superpartner produced at the LHC will 
cascade down to the long-lived /i NLSP. As the lightest Standard Model superpartner, the 
/i NLSP will then appear as a quasi-stable muon-like particle that can escape the collider 
detector before decaying into the gravitino; cf. [45, 40, 41, 46, 44] and references therein. 
Thus, one would find signatures that are very different from the excess in missing energy 
expected in the alternative neutralino LSP scenarios. 

The most dramatic catalytic enhancement is seen in the ^Li and ^Be production triggered 
by the formation of bound states of ^He with a (generic) negatively charged relic that 
we call X~. The catalytic path to ^Li and ^Be is shown by the following sequence of 



transformations [12, 31]: 

X- ^ (^HeX-) ^ ^Li, (1.1) 

X- -^ (^HeX-) ^ (^BeX-) ^ ^Be. (1.2) 

Although the (^HeX~) system has a binding energy of about 350 keV, its formation is 
delayed until T = 8 keV by an overwhelmingly large number of energetic photons that 
photo-dissociate (^HeX~). Thus, the (*^HeX~) bound state serves essentially as a "bottle- 
neck" for ^Li production, whereas the path to ^Be goes through the "double bottleneck" of 
(^HeX~) and (^BeX~). The key for the nuclear catalysis is an enormous enhancement of 
the reaction rates in the photonless recoil reactions mediated by X^ [12, 31]: 

(^HeX-) + D -^ ^Li + X- (1.3) 

{^BeX-)+n -^ ^Be + X^. (1.4) 

Indeed, since the rates of these catalyzed reactions exceed the SBBN rates for the production 
of ®Li and ^Be by many orders of magnitude, one finds a strong sensitivity of the efficiency 
of the primordial ^Li and ^Be production to the abundance of X~ at relevant times. For 
the case of X~ being a thermal relic, this abundance is governed by the X^ annihilation 
rate and by the X~ lifetime tx-- In turn, observationally inferred limits on the primordial 
abundances of both ^Li and ^Be will impose significant constraints on the lifetime of X~, 
its mass, and its interactions. The limits on the lifetime and the abundance are particularly 
interesting in view of the possible catalytic solution to the so-called ^Li problem [12, 21], 
which is a persistent discrepancy between the predicted primordial amount of ^Li and a 
factor of 2-3 lower observations of ^Li in the atmospheres of the metal-poor stars. 

For the gravitino LSP scenarios in which X~ is identified with a negatively charged 
slepton NLSP, X~ = /|f , the constraints imposed by primordial ^Li catalysis have already 
been analyzed in detail in a number of publications [12, 17, 18, 19, 20, 21, 22, 26, 27, 28, 
32, 33]. Assuming a standard cosmological history that leads to a typical thermal /i relic 
abundance, the bound from ^Li catalysis translates into an upper limit on the /i lifetime of 
Tj = Tx' ^ 5 X 10^ s [12, 20, 22, 26]. In collider-accessible regions of the parameter space, 
this limit is found to be considerably more restrictive than the BBN constraints associated 
with electromagnetic/hadronic energy release from /i decays [17, 18, 19, 22, 28, 32]. The 
Tj limit implies a gravitino mass mg well below 10% of the slepton NLSP mass mj for 
mj^ < 0{1 TeV) [18]. This seems to exclude a kinematical determination of mg [47] at 
the next generation of particle accelerators, which might have been feasible for O.lmj < 
iTiQ < mj [48, 49]. Consequently, the tj limit puts a big question mark over an idea of 
a collider test of supergravity via the microscopic determination of the Planck scale [47] , 
as well as the mg-dependent collider test of thermal leptogenesis [42] . Moreover, a post- 
infiationary reheating temperature above Tr ~ 10^ GeV and thereby even the viability of 
thermal leptogenesis with hierarchical right-handed heavy Majorana neutrinos seems to be 
disfavored by the ^Li constraint within the Constrained Minimal Supersymmetric Standard 
Model (CMSSM) for a standard cosmological history [19, 26, 28]. In the CMSSM, the 



splitting between rrir^ and m^ required to evade the r;?^ limit translates also into a lower 
limit on the gaugino mass parameter, which is assumed to take on a universal value 1711/2 at 
the scale of grand unification. Indeed, for the natural gravitino LSP mass range in gravity- 
mediated SUSY breaking scenarios, the cosmologically favored region can be associated 
with a mass range of the colored superparticles {e.g., a gluino mass of nig > 2.5 TeV) for 
which it will be very difficult to probe SUSY at the LHC [26, 28]. 

In this paper we analyze the constraints imposed on SUSY models by the catalysis of 
a primordial ^Be abundance. While — based on the results of Ref. [31] — we do not expect 
the constraints from ^Be to be considerably tighter than those from ^Li, this analysis is 
warranted for a number of reasons: 

1. Observations of ®Li are extremely difficult because its lines are not resolved spec- 
troscopically with respect to the lines of ^Li. The claim of a "^Li plateau" with 
metallicity [50] at ~ 10% of the '^Li abundance is being challenged in the recent pa- 
per [51], and a new re-analysis of ^Li data is warranted as some of the observations 
may turn out to provide only upper limits. In fact, even the value of an upper limit 
on the primoridal ^Li abundance is subject to discussions: Many papers adopt upper 
limits on primordial ^Li/H = neu/n-n within a range from 10^^^ to 10^^°. Unlike 
^Li, ^Be is firmly detected in a significant number of stars at low metallicity, and its 
observational status is not in doubt. For the latest data on the ^Be abundance in 
metal-poor stars, see, e.g., [52, 53, 54]. 

2. ^Li is more fragile than ^Li and would burn more efficiently at lower temperatures. 
Therefore, if there is a (yet unconfirmed) stellar mechanism (see e.g. [55]) that resolves 
the lithium problem, i.e., that depletes ^Li by a factor of two or three, ^Li would have 
been depleted by an even larger factor. Such a stellar mechanism, however, would 
affect ^Be less than either ^Li or ®Li since both ^Li and ^Li are more fragile than ^Be, 
which thereby provides a more robust bound on New Physics. 

3. The nuclear physics rates that enter in the calculation of ^Be catalysis are dominated 
by resonances. Given the wealth of experimental information on the ^Be resonances 
[56] , this may eventually allow for very reliable calculations of the catalytic rates. 

While stating strong bounds on primordial abundances/lifetimes of X^, a lingering 
question remains: Is there an "island" of allowed lifetimes around tx- ~ 10^ s? One 
point emphasized in Refs. [13, 15, 25] is that bound states of X~ with protons, {pX~), 
may have a significant impact on the primordial abundances of elements such as ^Li, ^Li, 
^Be (and ^Be). Another point is related to X~ decays with very energetic decay products 
and the associated possibility of an environment in which ^Li may be destroyed efficiently. 
The second point has already been addressed [17, 22, 32], and it has been found that 
late energy injection cannot suppress ^Li down to an acceptable level if it is significantly 
overproduced at 8 keV [17, 22, 32] Despite its importance, the issue of possible {pX~) 
catalysis remained largely unresolved. While the importance of charge exchange reactions 



Table 1: Properties of {pX^) and C^HcX^) bound states. For (^KcX^), the binding energy E\^ includes 
a finite charge radius correction [f2]. Bohr radii ae are quoted for ideaUzed Bohr-type bound states. 
The given recombination temperatures T,- are understood as the temperatures at which the corresponding 
photo-dissociation rate becomes equal to the Hubble rate. 



bound state E\^{keV) aB(fni) Ti.(keV) 

(^HeX-) -347 3.6 8 

{pX~) -25 29 0.6 



of (pX^) with ^He that may reduce the abundance of (pX^) was already mentioned in 
the early paper [13], subsequent publications either ignored this issue [15] or generally 
underestimated the impact of charge exchange reactions [25] . Recognizing its importance 
for the whole CBBN paradigm, we revisit the catalysis by {pX~) bound states. Indeed, 
in this paper, we put this issue to rest by reaching the conclusion that less than 10% of 
^Li and ^Be synthesized at 8 keV could possibly be affected by {pX~) catalysis, whereas 
typically suppression factors in excess of 100 are needed in order to evade the corresponding 
limits. This clarifies that the islands in the parameter region with large abundances/large 
lifetimes, which were suggested to remain viable in Ref. [25], cannot exist. 

This paper is organized as follows. In the next section we estimate the charge exchange 
reaction rates relevant for (pX~)-mediated catalysis and consider their consequences for 
lithium and beryllium. In Sect. 3 we analyze the bounds on the lifetime/abundance of X^ 
imposed by observations of ^Be in stellar atmospheres at low metallicities. The resulting 
constraints on SUSY models with the gravitino LSP and a charged slepton NLSP are given 
in Sect. 4. We reach our conclusions in Sect. 5. 



2 {pX ) catalysis and charge exchange reactions 

We begin the discussion of the possibility of (pX^)-mediated catalysis by reminding the 
reader of the basic properties of the {pX~) and ("^HeX") bound states. Table 1 lists the 
corresponding binding energies E\^, Bohr radii ob, and recombination temperatures T^, 
where Tj. is understood as the temperature below which the rate of photo-dissociation of a 
bound state becomes smaller than the Hubble expansion rate. As the table suggests, the 
recombination of {pX~) bound states becomes efficient only after the temperature drops 
below 1 keV which corresponds to an age of the Universe on the order of a few weeks or 
so, and we will assume in this section that the X~ lifetime tx- is large enough to allow for 
this recombination to happen. Clearly, the presence of a large number density of negatively 
charged particles nx- during the recombination with helium leads to an overproduction 
of ®Li and ^Be [12, 31] due to the nuclear catalysis at 8 keV. To be specific, we choose 



Yx- = nx- /n-Q = 10^^ (with ub denoting the baryon number density), tx- -^ oo, and 
using our previous results [12, 31, 26] determine the abundances of hthium and berylhum 
at T = 1 keV from CBBN production: 

6Li/H|T=ikcV ^ 8 X 10~« ; 9Be/H|T=ikcV ^ 3 x IQ-^" , (2.5) 

which is about three orders of magnitude above the observational bounds; cf. Sect. 3 below. 
Some of the synthesized ^Li and ^Be will be in bound states with X~ . 

Below T = 1 keV, the concentration of {pX~) is controlled by a rapidly diminishing 
photo-dissociation rate and by depletion through the charge exchange reaction [13]: 

(pX-) + ^He ^ (^HeX") + p. (2.6) 

This reaction may have a very large rate as its cross section is determined by the actual 
size of the {pX^) bound state that is of the order of ae — 30 fm (cf. Tab. 1). In fact, most 
of the recombined states {pX~) are immediately intercepted by the reaction (2.6) so that 
the resulting {pX~) abundance would remain quite small at all temperatures. 

From studies of charge exchange reactions of muons on hydrogen, it is known that the 
muon is captured into highly excited states that have large orbital momenta and large 
principle quantum numbers. The radii of these excited orbits of muonic hydrogen are 
comparable to the Bohr radius of ordinary hydrogen. In the case of (2.6), the capture 
would mainly proceed to the n = 3 and n = 4 levels of the (*^HeX~) bound state. To 
estimate the cross section for the reaction (2.6), we employ a semiclassical approximation 
in which the motion of helium is described by a classical trajectory while the proton is 
treated quantum mechanically. The large values of n and / of the resulting (^HeX~) bound 
states give some justification to this treatment. 

Calling R the separation between ^He and X^ (or, more generally, the separation be- 
tween X~ and the incoming nucleus of charge Z), we now investigate the R value at which 
the proton loses its ability to bind to X' . The one-dimensional slice of the proton potential 
energy in the field of X" and "^He, 

V(r) = -^ + ^ (2.7) 

r |r — K| 

is plotted in Fig. 1. The limit of f? ^ oo corresponds to an unperturbed binding of the 
proton to X" with a binding energy of E^, = —25 keV (cf. Tab. 1). For Z > 1 and finite R, 
the curve has a maximum at positive values of r referred to as Knax- As the ^He nucleus 
comes closer, R decreases. For R values below some critical distance Rd, the binding energy 
of the proton becomes positive so that the tunneling of the proton to r ^ +oo starts to 
become viable. For even smaller values of R, one can find another "critical" distance Rc2 at 
which the probability for the tunneling of the proton becomes comparable to 1 due to the 
fly-by of the ^He nucleus. In principle, this is not an easy quantum mechanical problem, 
and we use further simplifications to overcome that. In order to estimate Rd, we employ 




Figure 1: Potential energy of the proton in the field of X^ at r = and an incoming nucleus at r = —Rc2- 
The potential energy is plotted along the line connecting X^ with ''He (solid line) or ®Li (dashed line), 
respectively. As the distance between the incoming nucleus and X~ decreases, the potential well becomes 
more narrow, and the proton ground state energy level is pushed upward. The critical deconfinement 



distance R 



c2 



is defined as the distance at which the energy of the bound state found variationally 



using (2.8) becomes larger than the height of the barrier Vmax to the right of X . 

the variational calculation of the proton energy in the potential (2.7) by using the trial 
wave function for the ground state, 



?/^(/i, u) = exp[-(/i - u)R/{2a)] x (1 + uR/bY 



XS) 



where /j, and z/ are elliptic coordinates and a and b the minimization parameters. The 
coordinates are defined as /i = (ri + r2)/{2R) and z/ = (ri — r2)/{2R), where ri and r2 are 
the proton-nucleus and proton-X" distances, respectively. We calculate the energy of the 
ground state E™ as a function of the distance between X~ and the incoming nucleus R. We 
determine Rd from E™{Rci) = and estimate i?c2 from E™{Rc2) = Knax which describes 
the situation when even a metastable bound state simply cannot exist. The cross section for 
the charge exchange reaction is then approximated by the geometric one with the impact 
parameter p = Rc2, o = vri?^25 which essentially assumes a deconfinement probability of 1 
for R < Rc2- This approximation would fail for very fast incoming nuclei. For the case 
considered in this paper, however, the fly-by time is much longer than the period with which 
the proton orbits within (pX^). The results of our estimates are presented in Table 2. As 
can be seen from this table, a ^He-X~ distance of ~ 95 fm is sufficient to release the proton 
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Table 2: Deconfining distances and charge exchange reaction cross sections on the (pX ) target for 
incoming nuclei with different charges Z. The Rd values are given in units of the (pX~)-Bohr radius 

Cg and in units of fni. 



Z Rci at E = Rc2 at E = V„iax cr = nR'^2 ^^ bn 



1 


1.4 (40 fm) 


1.4 (40 fm) 


51 bn 


2 


3.7 (107 fm) 


3.3 (95 fm) 


280 bn 


3 


5.8 (167 fm) 


4.7 (135 fm) 


580 bn 


4 


7.9 (230 fm) 


5.7 (160 fm) 


850 bn 



from the bound state. Consequently, our estimate points to a very large cross section of 
almost 300 bn for the charge exchange reaction (2.6). For Z = 1 we have a cross section 
~ 27ra|, which compares well with the results for the charge exchange cross section in the 
case of muon-hydrogen scattering [57]. 

Using this cross section, we incorporate the charge exchange reactions in the network of 
Boltzmann equations and calculate the residual concentration of {pX~) in a wide tempera- 
ture range. In the limit of infinite lifetimes, tx- -^ oo, we find that the abundance of {pX') 
reaches its peak at around T = 0.7 keV. Its maximum abundance at these temperatures 
can be well approximated as 

in max / T /^ 

'''(r,X'\ . -7 Y: 



J^.4xl0-'(i^j. ,2.) 

where we made the safe assumption^ of Yx- ^ y^ne-, which ensures the linear scaling in 
(2.9). 

Reference [25] makes the somewhat surprising suggestion that even a tiny fraction of 
surviving {pX^) bound states may cause a significant reduction of the ^Li abundance, given 
the significant uncertainty in the nuclear rates. It is easy to see, however, that two types 
of processes are possible for the colliding ^Li-(pX~) system: the charge exchange reaction 
and the nuclear reaction, 

A: {j)X-) + %i -^ (^LiX~)+p, (2.10) 

B: {pX-) + %i -^ X-+^He + ^He. (2.11) 

At first, we completely ignore process A and concentrate on process B which destroys 
^Li. To illustrate our point, we will assign the maximal possible rate to process B, which 



^Focusing on a minimal MSSM particle content, Yx- is determined by the standard chemical decoupling 
of X^ from the primordial plasma for which Yx- ^ Yi^^c holds unless mx- ^ 0(4 TeV); see Sect. 4 for a 
detailed discussion of Yx- ■ 



is given by the unitarity bound in the s-channel, 



\^^/B = -J2-\" ) = / nn = „1 /9 > (2-12J 



m^Li m^uVrm^ Tg/^ 



where the last expression is calculated in units of A^j^^cm^s~^mol~^, and Tg = T/10^ K. 
Comparing the destruction rate to the Hubble rate at a fiducial temperature of Tg = 0.008, 
we find that the former is much smaller, 



/_-„,\max„max 
\au)^ l\pX-) 



H 



~0.03<1 (2.13) 

79=0.008 



for Yx- = 0.01. This tells us that at most 3% of the ^Li synthesized at 8 keV could 
potentially be affected by {pX~) states via process B, and therefore the whole issue of 
nuclear uncertainties is irrelevant given the strength of the charge exchange reactions. 

Even though we find that process B makes no impact on ^Li, it is still interesting to 
examine whether (2.12) corresponds to a realistic rate. As can be seen from Fig. 1, for 
reaction B to occur, the proton has to tunnel to the left through a distance of at least 
100 fm inside the Coulomb barrier of ^Li in order to trigger its decay to two helium nuclei. 
Notice that the field of X~ does not facilitate this tunneling in any way. Effectively, this is 
the same distance for tunneling that a free ~ 50 keV energy proton would have to overcome 
in ap + ^Li -^ ^He + ^He reaction. Therefore, one expects an exponential suppression of the 
corresponding probability by exp(— a/-Eg/50 keV) ^ 1, where Eq is the Gamow energy. 
At the same time, the deconfining rate is 100% as long as the impact parameter is less 
or equal to i?c2; leading to the inevitable conclusion that the rate of the charge exchange 
reaction A greatly exceeds that of the nuclear process B, 

^^"^^ « 1. (2.14) 



{av) 



Returning to the calculation of Ref. [25] that gives a very large estimate for process B, 
we believe that this estimate is probably an artefact of assuming a "frozen" profile for the 
proton wave function. In contrast, in the correct approach, the proton wave function is 
easily polarized and deconfined by the incoming heavier nucleus. Further doubts in the 
validity of the estimates in [25] are cast by the answer for the cross section at extremely 
small energies. For example, it follows from [25] that ctb has almost an atomic size cross 
section, whereas this is a nuclear reaction between objects of nuclear size. The key difference 
between our treatment of (pX~)-induced catalysis and that of Ref. [25] is the significant 
underestimation of the strength of the charge exchange reactions in the latter work, which 
in turn leads to overestimates in expressions analogous to our Eqs. (2.9) and (2.13). 

To conclude this section, we consider an interesting way of having an impact of {pX^) 
on ^Li. Indeed, a successive chain of charge exchange reactions can lead to molecular states 



that are finally destroyed in nuclear reactions with protons: 

{pX^)+^U^{'^UX~)+p (2.15) 

(pX-) + CUX-) ^ {%iX^) + p (2.16) 

(pX-) + CUX^) ^ (^LiXg-) +p (2.17) 

C'LiXg") +p-^^Be + 3X" or ^He + ^He + 3X~. (2.18) 

In the last step of this chain, (^LiXg") "ammonium" has a chance for a nuclear interaction 
with protons or helium unsuppressed by a residual Coulomb barrier since (^LiXg") is a very 
compact object. A similar chain exists for ^Be and ^Be where the sequence of the charge 
exchange reactions can proceed until Be-X~ "methane," {^BeX^). It is important to note 
that the efficiency of this chain reaction depends very sensitively on the concentration of 
the {pX^) bound states and on the mass of the X^ particle. The latter enters through 
the average relative velocity of two heavy objects, e.g., {pX~) and (^LiX~), which in turn 
scales as m-^l . Therefore, in the limit of an infinitely heavy X^ and with our treatment of 
the charge exchange reactions, the chain will be cut off right at the first step, terminating 
at (^LiX~). For weak scale relics, the suppression of the average velocity of X~-containing 
bound states relative to the velocity of light nuclei is from one to two orders of magnitude. 

Using the charge exchange rates estimated in this section, we run the set of Boltzmann 
equations to determine the residual concentrations of {pX~) and of the molecular bound 
states of ^Li with X~. The results are plotted in Fig. 2. As one can see, an initial concen- 
tration of X~ per nucleon of Yx- = 10~^ results in a (pX^) abundance that never exceeds 
the maximum (2.9), leading to a progressively diminishing number of molecular states. 

There is another plausible mechanism for the destruction of (^LiX~) that is related to 
the recoil of ^Li freed in the decay of X~ . Typical kinetic energies of ^Li after the decay 
are comparable to that on the orbit ~ 700 keV. In the center of mass frame with plasma 
protons, this corresponds to energies of ~ 100 keV. At such a center of mass energy, some of 
the ^Li nuclei released in the decays of the bound states will be destroyed. However, given 
that only a small fraction of ^Li is locked in bound states with X~and that the destruction 
rate is smaller than the thermalization rate, it is safe to conclude that also this mechanism 
cannot lead to large overall depletion factors for ^Li. The same argument applies to ^Be. 

To conclude this section, neither lithium nor beryllium synthesized in CBBN processes at 
8 keV would be affected in any significant way by the subsequent generation of (pX^) bound 
states. Thus, the part of the parameter space with a typical freeze-out X^ abundance and 
a long X~ lifetime is confidently ruled out, which is shown in more detail in the following 
sections. 



3 ^Be constraints on the X lifetime and abundance 

In order to constrain the {tx-jYx-) parameter space from the catalytic path (1.2) to 
^Be [31], we need to set an upper limit on its primordial abundance from existing observa- 
tions. It is generally accepted that the galactic evolution of the abundances of Be, along 



Yx- = 10-^ nix- = 100 GeV, r^- -^ oo 




0.001 



tg 



Figure 2: Evolution of primordial abundances as a function of time (or temperature Tg) from the input 
Yx- = 0.01, rnx- = 100 GeV, and tx- -^ oo. The solid (red in the web version) line of the (pX^) 
abundance reaches its maximum of ^ 4 x 10^^ at Tg ~ 0.008. Around the same temperatures, the yield 
of unbounded X~, Y^l'^, starts to decline more rapidly since it is removed by the recombination with p 
followed by the charge exchange reaction on ^He. 



with Li and B, are dominated by cosmic-ray nucleosynthesis. While Be is burned rapidly 
in stellar centers, it is produced in cosmic rays by the spallation reactions of fast protons 
and a particles hitting ambient CNO nuclei [58] . As a consequence, the abundances of Be 
and O are linked, leading to a secondary scaling, Be oc O^ [59]. On the other hand, inverse 
spallation reactions of CNO nuclei, both produced and accelerated in supernovae, will give 
a Be yield that is essentially independent of the metallicity of the interstellar medium. Such 
primary processes, leading to Be oc O, are expected to play a major role during the early 
galactic epochs [60]. 

The produced Be is subsequently supplemented in the outer layers of stars. Thus, 
old stars which are far from the galactic center (and thereby less affected by the galactic 
chemical evolution) bear the potential to encode any pre-galactic origin of Be. Indeed, 
Be has been observed in a number of Population II halo stars at very low metallicities 
[Fe/H] < —2.5; [A] = log^g^ + 12. Particulary noteworthy is the detection in the star 
G 64-12 at [Fe/H] ^ -3.3 [52]. The star's high Be value of logio(Be/H) ~ -13.05 might 
suggest a possible flattening in the Be trend during the early evolutionary phases of our 
galaxy [52]. Whether this really points to a primordial plateau or whether this indicates a 
Be dispersion at lowest metallicities [54] is not clear at present. 
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Figure 3a shows the original Be detection in the star G 64-12^ (filled dot) along with a 
subset of data points taken from Fig. 3a of Ref. [52] . The data of Fig. 3b are taken from 
Fig. 6b of the recent work [61] which also uses [0/H] as a metallicity tracer. The latter 
paper discusses the implications of a new temperature scale on the abundances of Li, Be, 
and B. In principle, different assumed physical parameters which characterize the stellar 
atmosphere may result in large systematic shifts of the inferred abundances. In this regard, 
it is important to note that Be is not overly sensitive to the assumed surface temperature 
of the halo dwarfs [61]. In the following we thus shall take a pragmatic approach: In 
both Fig. 3a and Fig. 3b, we obtain the least squares weighted mean (dashed lines) for a 
representative sample of stars at lowest metallicities. From the variance of the fit, we can 
extract a nominal 3o" upper limit (solid lines) on primordial ^Be. From Fig. 3b, we find 

logioBe/Hlhigh = -12.68 => ^Be/R < 2.1 x 10"^^ . (3.19) 

Conversely, Fig. 3a yields ^Be/H < 10^^^ while fitting only the last two data points with 
[0/H] < —1.3 in Fig. 3b would give ^Be/H < 1.3 x 10^^^. In our context, those values are 
less conservative so that we use (3.19) in the following. In Fig. 3b we have additionally 
fitted for a primordial component, ^Be/H|p, in combination with a primary scaling, ^Be/H = 
kO/H. It seems, however, that a purely primary mechanism with k ~ 2.9 fits the data 
best since ^Be/H|p comes out negligibly small. ^ Finally, we are aware that neither of the 
fitted mean values in Fig. 3 is very good in terms of x^- However, a firm conjecture of a 
Be plateau is not the purpose of this work, and indeed (3.19) does provide a sufficiently 
conservative limit to work with. 

We can now confront the constraint (3.19) with the CBBN yield of ^Be obtained by 
solving the associated Boltzmann equations. The central input parameter for the catalytic 
production of ^Be and ®Li is the abundance of X^ at the time of its recombination with 
■^He. Above 10 keV, we can track the resulting (^HeX~) abundance by using the Saha- 
equation since photo-dissociation proceeds rapidly. Only at T ^ 8 keV, (*^HeX~) starts 
to build up efficiently, and we couple it into the full set of Boltzmann equations. We use 
a ^He-X~ recombination cross section that is based on the work of Ref. [12]. It takes 
into account the finite size of the nucleus and includes a-captures into IS as well as 2S 
states. For the cross section of catalyzed ^Li production (1.3), we employ the result of a 
nuclear three-body calculation [20]. The path to ^Be proceeds via (^BeX^) bound states 
which are formed by the radiative fusion ^He -|- ("^HeX~) -^ (^BeX~) -|- 7. From there, 
^Be is subsequently produced by resonant neutron capture (*BeX~) + n —^ ^Be -|- X~ 
[31]; see Appendix A for the associated cross sections. In our code we can neglect the 
formation of (®BeX~) that proceeds via molecular bound states (^HeXg") [31]. This process 
becomes important only for a combination of large Yx- and large Tx-, i.e., a parameter 
region which is already excluded by ^Li overproduction. Also note that at the time when 



^For consistency with the rest of the data points, the ID LTE value has been plotted in Fig. 3a of the 
original reference [52] . 

^For a proper comparison between different assumed surface temperature scales and corresponding fits 
of primary versus secondary scaling, see Ref. [61]. 
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Figure 3: Observations of Be in Pop II halo stars. In the left panel (a), the data is taken from Fig. 3a 
of Ref. [52] and is plotted as a function of [Fe/H]. The right panel (b) shows the data from Fig. 6b of 
Ref. [61] where [O/H] provides the metallicity indicator. The filled dots depict the data points associated 
with the star G 64-12. The solid lines give the inferred nominal upper limits on ^Be from the weighted 
mean (dashed lines) of a sample of stars at lowest metallicity. Also shown in Fig. 3b is a fit of a primary 
scaling of Be; see main text. 



E, 



^BeX ) form, their photo-dissociation is not important because of the high binding energy 
('BeX-) ^ ^^gQ ^^^ ^^2]. Finally, for Tg < 0.2, the SBBN n abundance can already 
be tracked well by including the processes D + D — > n + ^He, T + D ^ n + ^He, and 
^He + n — > p + T in the reaction network [62]. Those cross-sections as well as the one from 
p- induced ^Li destruction can be found, e.g., in [63]. It is important to note that we assume 
the SBBN central value for the deuterium abundance. The early decays of X~ may result 
in an injection of nucleons into the system. This typically drives the deuterium abundance 
upward, resulting in an enhanced number of neutrons at later times and therefore in an 
increased output of ^Be, with the general scaling ^Be ~ const x (D/Dsbbn)^- We choose to 
disregard this effect, noting its model-dependent character. We are allowed to do so since 
its inclusion can only make the ^Be-derived bound on the X~ abundance stronger. 

We should mention at this point that anyone attempting a precision calculation of ^Be 
within the CBBN framework would have to include an additional channel related to the 
early production of beryllium as pointed out in [21]. This channel consists of X~ capture 
on ^Be with a subsequent p-induced reaction producing (^BX^), which then beta decays 
to the (®BeX~) bound state: 



X- 



^BeX- 



^BX- 



^BeX- 



^Be. 



(3.20) 



The efficiency of this chain is directly proportional to a rather small ^Be abundance. There- 
fore, the final output of ^Be/H via (3.20) is never very large, but could reach the level of 
~ 0(10^^^) for large abundances of X~~ . For the purpose of setting limits on particle physics 
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Figure 4: Evolution of catalyzed ^Li and ^Be production shown together with the formation of the 
"bottle-neck" abundances of CHeX") and (^BeX") for F^"^ = 5 x 10"'' and t^- = 5 x 10^ s. The dashed 
line gives the neutron abundance while the dotted line shows the abundance of free X~ . 

models, we are allowed to ignore this early chain (3.20), noting that it is generally subdom- 
inant and also mo del- dependent. In particular, the chain (3.20) depends on the properties 
of X~ [21], as well as on the non-thermal processes that can affect the ^Be abundance"^. 

Figure 4 shows the evolution of catalyzed ^Li and ^Be production from the solution of 
the corresponding set of Boltzmann equations below T = 10 keV. We parameterize Yx- 
by the X~ abundance prior to decay by introducing Y^^, where the superscript "dec" 
stands for decoupling, and by the X~ lifetime tx 
any moment during BBN is given by Yx- {t) 



the curves in Fig. 4, the values y^!f 



Y^°!r X exp 



so that the (total) X~ abundance at 
t/rx). In particular, to obtain 
5 X 10~* and tx- = 5 x 10'^ s are used. When 
the "bottle-neck" abundances of (^HeX~) and (^BeX") form, the catalytic paths (1.1) 
and (1.2) to ^Li and ^Be open up, resulting in the asymptotic values ^Be/H ~ 10~^^ and 
^Li/H :ii 3 X 10^^^. The dashed line shows the neutron abundance and the dotted line the 
free X~ abundance, which is dominated by its exponential decay. We remark in passing 



that residual recombinations of ^He with X lead to the crossing of the (^HeX 
lines at late time. 



and X. 



free 



In Fig. 5 we obtain exclusion boundaries from catalyzed ^Be and ^Li production in the 
{tx- iY^^) parameter space. For convenience of the reader, the X~ number density n'^1 



■^In addition, as pointed out in recent Ref. [75], the beta decay of ^B occurs predominantly to the excited 
states of '^Be, which will likely result in a break-up of (*BeX~) and further reduction of the efficiency of 
the chain (3.20). 
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Figure 5: Contour plot of CBBN abundance yields of ^Li and ^Be in the {tx- , Y^''^) plane. The sohd line 
shows the limit (3.19). The region above this line is excluded by ^Be overproduction. The lower (upper) 
boundary of the band corresponds to ^Li/H = 10"^^ (IQ-^"). The y-axis on the right-hand side indicates 
the X~ number density n'^'i normalized to the entropy density s. The cross shows the parameter point 
considered in Fig. 4. 

normalized to the entropy density s is given on the j/-axis on the right-hand side. Above 
the sohd hne, ^Be is in excess with respect to (3.19) and thus excluded. The shown band 
reflects the uncertainties in the observational determination of ^Li. On the lower border, 
^Li/H = IQ-ii is fulfilled while ^Li/H = 10-i° holds on the upper border of the band. The 
cross indicates the exemplary parameter point considered in Fig. 4. At large lifetimes, the 
linear scaling of ^Li with Yx- can easily be seen from the boundaries of the band. Note 
that we find ^Be/®Li in the interval between IQ-^ and IQ-^, whenever CBBN is efficient, 
which confirms the observation already made in Ref. [31]. 



4 Implications for Supersymmetric Models 

Let us now address the implications of the results derived above to SUSY extensions of 
the Standard Model in which the gravitino G is the LSP and a charged slepton Zi the 
NLSP [37, 34, 39, 40, 35, 36, 41, 17, 18, 19, 26, 27, 28, 33]. As the spin-3/2 superpartner 
of the graviton, the gravitino is an extremely weakly interacting particle with supergravity 
couplings [64, 65] that are suppressed by the (reduced) Planck scale [66] Mp = 2.4 x 
IQi^ GeV. Thereby, the negatively charged l^ can be the long-lived X~ with a lifetime 
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Figure 6: Cosmological constraints on the masses of the gravitino LSP and a charged slepton NLSP for 
a yield Yf'''' given by (4.22). The gray band indicates n^^€ n^y^. Above this band, 17g > 0.126. On the 

thin sohd hne labeled with / values only / floM is provided by fi -"""^ . The dotted lines show contours of 
Tj . Due to CBBN, the region below the solid and the long-dash-dotted (red in the web version) lines is 
disfavored by observationally inferred abundances of ^Be and ^Li, respectively. The effect of electromagnetic 
and hadronic energy injection on primordial D disfavors the regions inside the short-dash-dotted (blue in 
the web version) curves and to the right or inside of the short-dashed (blue in the web version) curves, 
respectively. The region below the dashed (green in the web version) line is disfavored by the effect of 
electromagnetic energy injection on ^He/D. While the constraints from hadronic energy injection are 
obtained for a purely 'right-handed' li ~ /r NLSP, the ones from electromagnetic energy injection are valid 
for the Ti NLSP case with a visible electromagnetic energy of -Evis = Ecm = O.S-Et- released in ti -^ Gt. 



Tx = Tj governed by the decay Zi 
rr ^ T~\h - 



Gt,^ 



Gl) 



AMm^Ml 



G 



m- 




(4.21) 



Indeed, tj > 10^ s occurs in a large region of natural values of the gravitino mass mg and 
the slepton mass mj , as illustrated by the tj -contours (dotted lines) in Fig. 6. 

For a standard cosmological history with a post-inflationary reheating temperature Tr 
above the decoupling temperature of the /i NLSP, Tf < mj /20 [37], the /i NLSP freezes 



out of the primordial plasma as a cold thermal relic so that its yield after decoupling Y~ 

^1 



dec ; 



IS 



^We assume R-parity conservation. For the case of broken R-parity, see e.g. [67]. 
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governed by its mass and its annihilation rate. Thereby, 1^^'^'^ becomes sensitive to the mass 
spectrum and the couphngs of the SUSY modeL In this work, we work with a representative 
yield that is quite typical for an electrically charged massive thermal relic [37, 39, 19]^ 

ydcc ^ _«! = 2 Y^^ = 0.8 X 10"^ I j^ — ) , (4.22) 

'i ne 'i VlOOGeV/ ^ ^ 

where nj denotes the total /i number density assuming an equal number density of pos- 
itively and negatively charged /I's. Note that the yield (4.22) is in good agreement with 
the curve in Fig. 1 of Ref. [37] that has been derived for the case of a purely 'right-handed' 
Ti ~ tr NLSP with a mass that is significantly below the masses of the lighter selectron 
and the lighter smuon, m^^ ^ me^jxi, and with a bino-like lightest neutralino, Xi — B, 
that has a mass of m^ = 1.1 m^-^. In the case of an approximate slepton mass degeneracy, 
"^n ^ ^TT'ei,jii ^ l.lm^-^, the Ti NLSP yield (4.22) can become twice as large due to slepton 
coannihilation processes [37, 19]. Approaching the Xi^^i coannihilation region, m^o ~ mf^, 
even larger enhancement factors occur; see e.g. Fig. 3 in Ref. [19]. On the other hand, 
a sizable left-right mixing of the stau NLSP is associated with an increase of its MSSM 
couplings and thus with a reduction of Yj . Moreover, an exceptional reduction of Yj can 
occur also in a non-standard thermal history with late-time entropy production after the 
decoupling of the /i NLSP and before BBN [69, 19, 20] or in low Tr scenarios [23]. Noting 
that both cases require substantial modifications to the MSSM field content at or below 
the weak scale, we disregard such possibilities and focus in the remainder of this work on 
the more generic Yj values described by (4.22). 

Confronting (4.22) with our limits shown in Fig. 5, we obtain the CBBN constraints 
shown in Figs. 6, 7, and 8. In each figure, it is the region to the right of the long-dash-dotted 
(red in the web version) lines and the one to the right of the solid line in between those lines 
that is disfavored by ^Li/H|p < 10~^^ and 10"^*^ and by the ^Be limit (3.19), respectively. 
While the emphasis is on the limits from CBBN of ^Be (and ^Li), the following additional 
cosmological constraints are shown for comparison: 

• Since each /i NLSP decays into one G LSP, these decays lead to a non-thermally 
produced (NTP) gravitino density [37, 40]: 

Qf^h' = mg Ff ^ nB(To)/iVpc , (4.23) 

where pc/[^B(To)/i^] = 42.1 GeV [66]. This contributes to the relic gravitino density 
Qq which should not exceed the observationally inferred dark matter density ^dm and 
thus imposes an additional constraint on the model. Accordingly, we show shaded 
regions in Figs. 6, 7, and 8, in which the Q%^^h'^ values obtained with (4.22) agree 
with the nominal 3a range of ^dmh'^ inferred with a restrictive six-parameter "vanilla" 
model from the three year data set of the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite [1] 

nl^^h' = 0.1051^30 (4.24) 



^For a recent thorough study of the decouphng yield of a charged rehc, see Ref. 
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with h = 0.73^003 denoting the Hubble constant in units of lOOkmMpc^^s^^. In each 
figure, the parameter space above the shaded region is disfavored by the dark matter 
constraint Q%^^ < Qbm- With any additional contribution to ^dm — such as an axion 
density or a thermally produced gravitino density Qjr^ — this constraint can become 
even more restrictive. In Fig. 6 this is indicated by the thin solid contours labeled 
with / = 0.1, 0.01, 10-^ and 10-^ on which (4.23) obtained with (4.22) satisfies 
/ngTP ^ 0.126, respectively. 

In a /i NLSP decay. Standard Model particles are emitted in addition to the gravitino 
which can affect the abundances of the primordial light elements. While the associ- 
ated hadronic/electromagnetic energy release seems to affect the CBBN constraints 
from ^Li and ^Be only mildly — as discussed in the Introduction — it could alter in a 
substantial way the primordial fractions of D/H and ^He/D [40, 35, 41, 17, 22]. The 
effect of hadronic energy injection on primordial D disfavors the regions inside the 
short-dash-dotted (blue in the web version) curves shown in Figs. 6, 7, and 8.^ These 
curves are obtained from the upper limits on Y^'^'^ that are given in Fig. 11 of Ref. [41] 

as derived from a computation of the 4-body decay of a purely 'right-handed' /i ~ /r 
NLSP into the gravitino, the tau, and a quark-antiquark pair [41]. Note that these 
upper limits are based on the severe and conservative upper bounds on the released 
hadronic energy (95% CL) obtained in [9] for observationally inferred values of the 
primordial D abundance (see references cited in [9]): 

(D/H)mcan = {2.78to.tt) X 10"^ ^ severe constraint, (4.25) 

(D/H)high = (3.98+°j^) X 10"^ ^ conservative constraint. (4.26) 

Without trying to give extra credence to a rather high value of D/H in (4.26), follow- 
ing [9], we simply take it as a limiting value for D/H. 

The regions disfavored by electromagnetic energy injection are shown in Fig. 6 only. 
Here it is the region to the right or inside of the short-dashed (blue in the web version) 
curves and the region to the right of the long-dashed (green in the web version) line 
that are disfavored by the primordial abundances of D and ^He/D, respectively. These 
curves are obtained for the stau NLSP case h = ti, i.e., for a 'visible' electromagnetic 
energy of -Evis = Cem = 0.3 Er of the tau energy Er = (m?^ — m~ + m1)/2mf^ released 

in Ti -^ Gt,^ where the D^^ and ^He/D constraints result from the Inlsp limits given 
in Fig. 42 of Ref. [9] and the D^°°^ constraint from the Vnlsp limit given in Fig. 6 of 



''Additional constraints on hadronic energy release are imposed by the primordial abundances of *He, 
■^Hc/D, ^Li, and ^Li/^Li [5, 6, 8, 9, 10, 17]. However, in the region allowed by the ^Be and ^Li constraints 
from bound-state effects, i.e., Tf^ < 10'* s, the considered D constraint on hadronic energy release is the 
dominant one as can be seen e.g. in Figs. 38-41 of Ref. [9] and in Figs. 6-8 of Ref. [10]. 

*For the selectron NLSP, li = ei, with E'vis = eom = {m~ — m~ + rn^)/2m'^^ given by the 'full' electron 

energy released in ei -^ Ge, the bounds become more severe as shown explicitly in Fig. 12 of Ref [41]. 
However, a comparison with Fig. 5 shown above shows that the constraints from electromagnetic energy 
injection will still be significantly less restrictive than the CBBN constraints from ^Be and ^Li. 
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Ref. [7]; see also Fig. 9 (lower panel) in Ref. [41]. As noted before, the elevated content 
of D leads to the enhancement of CBBN-produced ^Li and ^Be. For example, if non- 
thermal processes boost the deuterium abundance to the level of (4.26), it would lead 
to an enhancement of the ^Li output by a factor of ~ 2, while the corresponding 
enhancement factor in the case of ^Be is about 4. 

Simple comparison shows that the ^Be constraint (together with the one from ®Li) 
provides the most restrictive upper limit on m^ for a given mj in the collider-accessible 
region below 1 TeV. Indeed, in the most conservative case, / = 1, the dark matter con- 
straint Qj~^^ < Qdm disfavors {m^,mfj combinations associated with mg > 200 GeV and 
^ ^ 800 GeV, which is a mass range that will be difficult to probe at the LHC or at 
the International Linear Collider (ILC). Moreover, the 'electromagnetic' Dcm and ^He/D 
constraints are always less restrictive than the CBBN constraints from ^Be and ^Li. Only 
the 'hadronic' constraint Dhad can potentially compete with the CBBN constraints. For 
instance, this occurs for (mg, -my- ) combinations associated with mg > 60 GeV and a heavy 

/i NLSP with mj > 900 GeV. Accordingly, the presented CBBN constraints from ^Be and 
^Li are the most relevant ones in the mj range that will be accessible at the next generation 
of particle accelerators. 

For mr below 1 TeV, the new ^Be constraint can also be considered as the most robust 
BBN constraint. Indeed, the difference between the contours labeled with D'^™^'^ and D*^*^^ 
demonstrates that the D constraints are associated with a significant uncertainty related 
to the assumed upper limit on the primordial D/H fraction; cf. (4.25) and (4.26). The 
uncertainty associated with assumed upper limits on the primordial ^Li/H fraction, which 
can differ by (even more than) an order of magnitude, is indicated by the difference between 
the ^Li contours obtained for ^Li/H|p < 10^^^ and 10^^°. Indeed, since the observational 
status of ^Be is in better shape than the one of ^Li and since ^Be is less fragile than ^Li and 
^Li and thus less affected by stellar mechanisms, the ^Be constraint (which is represented 
by a single line) can be considered to be more robust than the ^Li constraint. 

Here we would like to emphasize that the ^Be and ^Li constraints are the ones that are 
the least sensitive to the precise value of Y^^'^ in the region Y^'^'^ > 10^'^. This results from 

the fact that the limits are very steep in that region, as can be seen in Fig. 5. Indeed, a 
yield that is twice as large as (4.22) will affect the position of the ^Be and ^Li constraints 
only very mildly. In contrast, such an enhanced yield — as encountered, e.g., in the case 
of slepton coannihilations — leads to significant changes of the dark matter constraint and 
the BBN constraints associated with hadronic/electromagnetic energy injection, as can be 
seen explicitly in Fig. 17 of Ref. [18]. It should also be noted that an elevated slepton 
yield can lead to an additional non-thermal output of ^Li for tj > few x 10^ s. This 
is because energetic spallation debris of destroyed ^He nuclei from slepton decays can hit 
ambient *^He and thereby fuse ®Li [8, 9, 10]. This mechanism depends sensitively on the 
hadronic branching ratio B^ of the 4-body slepton decay into the gravitino, the associated 
lepton, and a quark- ant iquark pair for which typically -Bh ^ 3 x 10^'^ for mj < 2 TeV 
(see Fig. 5 of Ref. [41]). Indeed, as discussed in Ref. [29], for those branching ratios, the 
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Figure 7: Contours of mg (dotted lines) as a function of tj- and mj . Assuming a slepton yield Y-'^'^ 
given by (4.22), constraints from CBBN of ^Be and ^Li are obtained as shown by the solid line and the 
long-dash-dotted (red in the web version) lines, respectively. For a purely 'right-handed' li ~ /r NLSP with 
a yield (4.22), BBN constraints from effects of hadronic energy injection on D are obtained as indicated by 
the short-dash-dotted (blue in the web version) lines. In the shaded region, Q^^h? agrees with (4.24). 



effect of CBBN on ^Li is the dominant one in the region which is not aheady excluded by 
the D constraint. Thus, for mj ^1.5 TeV, our obtained hmits on ^Li overproduction are 

only marginally affected by the hadronic energy release of /i-decays. However, for larger 
slepton masses, i.e., for scenarios of large Y-^'^ in conjunction with B\^ > 10~^, the hadronic 

production of ^Li becomes efficient so that only a simultaneous treatment of both effects 
can decide on the accurate ^Li BBN output. Note that this can make our presented limits 
on ^Li only stronger. Thus, we are on the conservative side when neglecting such additional 
contributions. 

Considering the CBBN constraints in Figs. 6, 7, and 8, one finds that the constraints 
from ^Be and ^Li are in close proximity of each other. This coincidence is of course related 
to the fact that the ratio of maximally allowed values for ^Be and ^Li is in rough agreement 
with the efficiency of producing ^Be and ®Li per each long-lived negatively charged slepton. 
Given the close proximity of these limits, it is tempting to speculate about a possible CBBN 
origin of primordial abundances of ^Li and ^Be at lowest metallicities [31]. This possibility 
applies to the gravitino dark matter scenarios that are located at or only slightly to the left 
of these constraints. It will be interesting to see whether also a solution of the ^Li problem 
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Figure 8: Contours of T^^"" imposed by fl^^h^ + Q^^^h^ < floMh^ < 0.126 for c = 1 (dashed lines) and 
c = 7 (dotted lines) as a function of tj and mj . The other cosniological constraints are identical to the 
ones shown in Fig. 7. 

can be found in these scenarios. Indeed, several ways to solve the ^Li problem have been 
proposed that could be relevant in this region [8, 36, 12, 17, 21, 25, 30]. However, a definitive 
answer will require an elaborate treatment of BBN in which all relevant effects from bound- 
state formation and electromagnetic/hadronic energy release are included simultaneously. 

Having discussed the generic features of the ^Be constraint and its comparison with other 
BBN constraints, we now would like to address its implications. Our task is facilitated by 
the ^Be constraint being very close to the one from ®Li, which makes those implications 
similar to the ones of the ^Li constraint [12, 17, 18, 19, 20, 22, 26, 27, 28, 33]. In this 
respect. Sect. 2 becomes important in which we show that the possibility of allowed islands 
in the parameter region with large Yj /large tj — which was advocated to remain viable in 
Ref. [25] — does not exist. Our present work does thereby reassure the conclusions drawn 
from the ^Li constraint in a decisive way : 



1. The gravitino mass m^ is constrained to values well below 10% of the slepton NLSP 
mass mj for mj < (9(1 TeV). This can be read off conveniently from the CBBN 
constraints shown in Fig. 6 and also from Fig. 7 in which rriQ contours (dotted lines) 



are given as a function of tj 



and mr . 



In particular. 



0.1 mr 



< 



mg < 



mr 



and 



thereby the kinematical determination of m^ proposed in [47] remains cosmologically 
disfavored at least at the next generation of particle accelerators. 
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2. The CBBN constraints disappear for a gravitino mass of m^ < 200 MeV provided 
mj > 80 GeV as supported by the non-observation of long-hved charged sleptons 
at the Large Electron Positron Colhder (LEP) [66]. This can be seen in Fig. 6. 
Accordingly, for gauge-mediated SUSY breaking leading to small values of rriQ, the 
CBBN constraint can be irrelevant. However, for mg > 10 GeV, as obtained in 
gravity- mediated SUSY breaking, the CBBN constraints impose a lower limit of mj > 
400 GeV as can be seen in Figs. 6 and 7. With CMSSM relations between the 
masses of the superparticles, this translates into a lower limit on the gluino mass 
TTig > 2.5 TeV [26, 27, 28]. Thereby, the ^Be constraint also points within the CMSSM 
to a cosmologically favored mass range that will be difficult to probe at the LHC. 

3. The ^Be constraint imposes an upper limit on the lifetime tj that ranges between 
6 X 10^ s at rriY = 100 GeV and 3 x 10^ s at mj = 1.5 TeV for a charged slepton 
NLSP with Y^^'^ described by (4.22), as can be seen in Fig. 7. Indeed, when the value 
of Y^'^'^ increases by a factor of 15, the tj limit decreases by only a factor of two. 

This mild mj dependence of the tj limit reflects the fact that the ^Be constraint 
is quite insensitive to the precise value of Yj . While the mg contours and the D 
constraints shown in Fig. 7 are specific to the gravitino LSP scenario with unbroken 
R-parity, the CBBN limits on tj shown in Fig. 7 apply to any scenario with a long- 
lived charged slepton described by the yield (4.22) and can thus be relevant for axino 
LSP scenarios [70, 71, 44] and scenarios with R-parity violation [67] as well. 

4. Gravitino dark matter can originate not only from NLSP decays but also from thermal 
scattering of particles in the hot primordial plasma. Thereby, the relic gravitino 
density Qg receives an additional contribution Qj? that depends basically linearly on 
the reheating temperature Tr after inflation [38, 41, 42, 43]. In turn, the dark matter 
constraint, Qj^ + Q^^^ < J^dm, can be translated into a conservative upper limit [33] 



2.37 X 10^ GeV f ^puh' - l ^g^^/i^ 

C2 I 0.1 



ra< ^'^'\T I .,° ]{^n,^i^)-Tr-. («7) 




which depends on the ratio of the gluino mass m-g and the li NLSP mass, c = 
m-g/mj > 1, at the weak scale. In Fig. 8, contours of T^'^^ obtained with Q^uh'^ ^ 
0.126 and Q%^^ given by the yield (4.22) are shown for the limiting case c = 1 (dashed 
lines) and for the case c = 7 (dotted lines), which is typical for universal soft SUSY 
breaking parameters at the scale of grand unification. Thus, also with the ^Be con- 
straint and with f2~'^^ included, c = 1 (c = 7) is found to be associated with a 
2^max yg^i^g of 3 X 10^ GeV (10^ GeV) in the cosmologically favored region. Indeed, 
the Tr constraints for c = 7 are consistent with earlier findings within the CMSSM 
and other constrained scenarios [19, 26, 72, 27, 28]. Since the ^Be constraint imposes 
Tj < 6 X 10^ s as discussed above, associated upper limits on the mass ratio c for 
a given lower limit on the reheating temperature — such as Tr > 10^ GeV required 
by successful thermal leptogenesis with hierarchical heavy right-handed Majorana 
neutrinos — can be inferred from Fig. 3 of Ref . [33] . 
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Let us conclude this section by considering prospects at future colliders. If the gravitino 
LSP scenario with a not too heavy charged slepton NLSP is realized in nature, the produc- 
tion and analysis of the (quasi-) stable charged sleptons will be a realistic option. Thereby, 
collider measurements oimj will become available [45, 46]. Moreover, with an experimental 
reconstruction of at least some part of the SUSY model, one will be able to calculate Y^'^^ 
reliably for a standard thermal history with Tr above the decoupling temperature Tf of the 
slepton NLSP. Confronting the obtained 1-^*^'^ with the CBBN constraints from ^Be and ^Li 

shown in Fig. 5 can then provide an upper limit on the lifetime tj of the long-lived slepton. 
Assuming the gravitino LSP scenario, this tj constraint together with the measured mj 
will imply an upper limit on the gravitino mass mg as can be seen in Figs. 6 and 7. In 
addition, if tj can be measured, e.g., by analyzing /i decays in a collider detector [45, 48] or 
in some additional stopper material [73, 74, 71, 49], one should be able to compare the ex- 
perimentally determined combination {tj , mj ) with the CBBN constraints shown in Fig. 7 
(without any assumption on the gravitino LSP scenario). While a finding of {tj ,mj) in 
the region disfavored by CBBN could point to a non-standard cosmological history with 
late-time entropy production [69, 19, 20] or to a low reheating temperature [23], it would 
be most remarkable to find {tj , mj ) in close vicinity of the CBBN constraints. Notwith- 
standing a rather large number of "if" s compounded in the previous sentences, the collider 
measurements could provide in this way an independent test for a hypothesis of the CBBN 
origin of ^Li and ^Be at lowest metallicities. Complementary to that it will be exciting to 
see new analyses of ^Be and ^Li data from future astrophysical observations. 



5 Conclusions 

From observations of beryllium in Population II halo stars at very low metallicities, we have 
extracted a nominal upper limit on primordial beryllium of ^Be/H < 2.1 x 10~^^. This limit 
allows one to set interesting constraints on models in which the primordial A = 8 divide 
is bridged by catalytic effects. Considering the primordial catalysis of ^Be via bound-state 
effects of a negatively charged massive relic X~ [31], we have derived rx--dependent upper 
limits on the X~ yield prior to decay, Y^'^f. For a typical relic abundance Y^^J^ ^ 3 x 10~^ 
(10~^), we find that this ^Be limit translates into an upper limit on the X~ lifetime of 
Tx- ^ 6 X 10^ s (10*^ s), which is quite comparable with the tx- limit inferred from the 
primordial catalysis of ^Li. Moreover, in the region where CBBN is efficient, we confirm 
that the ratio of the synthesized elements of ^Be/^Li lies in the range 10~^ — 10~^ [31], which 
provides perhaps the most model-independent prediction in the whole CBBN paradigm. 

We have clarified that the presence of {pX~) bound states cannot relax the Y^'^f limits 
at long lifetimes tx- in any substantial way. Indeed, we have shown explicitly that late-time 
effects of {pX~) bound states can affect the lithium and beryllium abundances synthesized 
at T ~ 8 keV by not more than 10%. Any substantial formation of {pX~) at T ~ 0.7 keV 
is immediately intercepted by the very efficient charge exchange reaction of (pX^) with 
*^He. This comes as no surprise given the large size of the {pX~) system ~ 30 fm and 
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the fact that the proton deconfinement probabihty approaches unity already for a ^He-X^ 
distance of ~ 95 fm. In particular, we find that the fractional density of protons in bound 
states does not exceed the level of ~ 10"^ for Yx- ^ ^4He- Correspondingly, even with a 
(pX~)-induced ^Li-destruction-cross section as large as the unitarity limit, at most a few 
percent of the synthesized ^Li could be destroyed. By the same argument, the ^Be yield 
also remains unaffected by late-time catalysis. Thus, we find that the possibility of allowed 
islands in the parameter region with typical Y^^J^ and large Tx- — which was advocated in 
Ref. [25] — does not exist. 

Applying the rx--dependent upper limits on Y^'^f derived from the primordial catalysis 
of ^Be, we have analyzed the new ^Be constraint in SUSY models in which the gravitino 
is the LSP and a long-lived charged slepton the NLSP, /i = X~ . For typical values of the 
slepton NLSP yield after decoupling, the ^Be constraint obtained in this paper is found 
in close vicinity to the constraint from the primordial catalysis of ^Li. Accordingly, the 
implications of the ^Be constraint for SUSY models do not differ much from the case of ^Li. 
The important advantage of ^Be-derived constraints is due to the fact that ^Be is firmly 
detected in a significant number of stars at low metallicities, while the status of ^Li obser- 
vations remains somewhat questionable. Another great virtue of ^Be-derived constraints is 
due to the expectations that an ^Be abundance would remain less affected compared to the 
one of ^Li by any stellar mechanism that might have caused the depletion in ^Li. Therefore, 
one would not expect any serious depletion factors between the (hypothetical) primordial 
fraction of ^Be and the actual observationally determined abundances of ^Be. 

One could question the calculational status of the CBBN chain (1.2) leading to ^Be. 
The rates for the first two reactions, resulting in (^HeX~) and (^BeX~), are determined by 
electromagnetic interactions and thus are not associated with large nuclear uncertainties. 
The final step on the way to ^Be is the neutron capture by the (®BeX~) bound state. It is 
dominated by a resonant transition, and the associated rate is believed to have a factor of 
a few uncertainty [31]. However, given the wealth of existing experimental data on the ^Be 
resonances, and rapid progress in nuclear calculations of few nucleon systems, one could 
hope that reasonably precise calculations of the catalytic rates may become available. 

The main shortcoming of our analysis is that only the catalytic effects below 10 keV 
are taken into account. In a generic framework of a hypothetical X' particle, this is fully 
justified. However, in the specific SUSY model with gravitino LSP/slepton NLSP, one could 
go one step further and by combining the catalytic effects with the energy injection effects. 
Fortunately for us, the question of ^Be synthesis is somewhat decoupled from the question 
of energy injection. The connection is mainly due to the modification of the deuterium 
abundance: For example, hadronic energy input at T ~ 30 keV would lead to a larger 
deuterium abundance, which in turn would enhance the neutron abundance resulting in 
a larger abundance of ^Be. Thus, the energy injection would lead to abundances of ^Be 
that are somewhat larger than the ones that we determined by using the standard input 
for the deuterium abundance. Since the main idea of our paper is to derive a conservative 
limit on SUSY models from ^Be, we can stay on the conservative side and disregard this 
modification, which would make the ^Be limit only stronger. 
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In comparison to other cosmological constraints on gravitino LSP scenarios with a long- 
hved charged slepton NLSP, we find that the ^Be constraint (together with the one from 
^Li) is the most relevant one in the colhder-accessible region of slepton masses below 1 TeV. 
Indeed, if a SUSY scenario with a long-lived /i is realized in nature, one might be able to 
determine the combination {tj , mj ) at collider experiments. Assuming standard cosmolog- 
ical history, one might also "invert" the collider data and infer Y^'^'^. Independently of the 
assumption of the gravitino LSP, these quantities can then be confronted with constraints 
on the {tj , Y^'^^) parameter space imposed by the primordial catalysis of ^Li and ^Be. 
Thereby, the CBBN constraints can be considered as predictions that could be tested in 
upcoming high-energy experiments. It will be most remarkable if collider measurements 
point to a {tj , mj ) combination in the vicinity of the CBBN constraints. Indeed, this 
could provide an experimental hint for the primordial catalysis being the origin of existing 
abundances of both ^Be and ^Li at lowest metallicities. However tenuous the BBN-LHC 
connection may seem at the moment, we expect a lot more clarity brought to this issue in 
the coming years. 

Note added - After the submission of this paper, a dedicated nuclear physics study 
of some CBBN reactions has appeared, Ref. [75]. It supports the conclusion of this paper 
about the large rate for the charge exchange reactions that remove {pX~). At the same time, 
this work finds non- negligible shifts, 0(100 keV), of the resonant energy levels employed 
in the ^Be production chain. This may affect the overall efficiency of ^Be production, and 
further investigations of the nuclear rates are needed. 
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A CBBN reaction rates below 10 keV 

In the following, we collect the key reaction rates, Nfi^{av), used in the numerical solutions 
of the Boltzmann equations. They are given in units of cm^s~^mol~^ and Tg = T/10^ K. 
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Recombination and photo-dissociation of X : 

^He + X- ^ (^HeX-) + 7 : 7900 Tg^^^ 

(^HeX-) + 7bg ^ ^He + X" : 1.85 x 10^° Tg^ exp (-4.O3/T9) 

p + X- ^ {pX~) + -f : 3980 Tg"^/^ 

{pX') + 7bg -^ p + X~ : 1.18 X 10^ T9-2 exp (-O.29/T9) 

Charge exchange reactions: 

(pX-) + ^He -^ (^HeX") + p : 3.9 x 10^° Tg^^ 

(pX-) + ^Li -^ (^LiX-) + p : 6.45 x 10^° Tg^^ 

(pX") + (^LiX-) ^ (6LiX2") +p : 3.37 x 10^ Tg^^ (1 TeV/mx-)^/^ 

(pX-) + (^LiX2-) -^ CUX^) + p : 5.25 x 10® Tg^^ (1 TeV/mx-)^/^ 

^Li and ^Be catalysis (from [20] and [31], respectively): 

(^HeX") + D^^Li + X" : 2.37 x 10® (1 - 0.34 Tg) Tg~^/^ exp (-5.33 Tg^^^^) 

^He + (^HeX") -^ (^BeX") + 7 : 10^ T'^^^ [ 0.95 exp (-I.O2/T9) 

+ 0.66 exp (-I.32/T9) ] 

(®BeX-)+n^^Be + X- : 2x10^ 



References 

[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 
[arXiv:astro-ph/0302209] . 

D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007) 
[arXiv:astro-ph/0603449] . 

E. Komatsu et al. [WMAP Collaboration], arXiv:0803.0547 [astro-ph]. 

[2] For a review of BBN constraints, see e.g.: S. Sarkar, Rept. Prog. Phys. 59, 1493 (1996) 
[arXiv:hep-ph/9602260] . 

[3] D. Lindley, Astrophys. J. 294, 1 (1985). 

J. R. Ellis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259, 175 (1985). 

R. J. Scherrer and M. S. Turner, Phys. Rev. D 33, 1585 (1986) [Erratum-ibid. D 34, 
3263 (1986)]. 



25 



[4] M. H. Reno and D. Seckel, Phys. Rev. D 37, 3441 (1988). 

S. Dimopoulos, R. Esmailzadeh, L. J. Hall and G. D. Starkman, Nucl. Phys. B 311, 
699 (1989). 

Yu. L. Levitan, I. M. Sobol, M. Y. Khlopov and V. M. Cliechetkin, Sov. J. Nucl. Phys. 
47, 109 (1988) [Yad. Fiz. 47, 168 (1988)]. 

[5] G. Sigl, K. Jedamzik, D. N. Schramm and V. S. Berezinsky, Phys. Rev. D 52, 6682 
(1995) [arXiv:astro-ph/9503094]. 



[6 

[7; 

[8 
[9 

[10 
ill 



[12 
[13 

[14 

[15 
[16 

[17; 

[18 
[19 



K. Jedamzik, Phys. Rev. Lett. 84, 3248 (2000) [arXiv:astro-ph/9909445]. 

R. H. Cyburt, J. R. Ellis, B. D. Fields and K. A. Olive, Phys. Rev. D 67, 103521 
(2003) [arXiv:astro-ph/0211258]. 

K. Jedamzik, Phys. Rev. D 70, 063524 (2004) [arXiv:astro-ph/0402344]. 

M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 71, 083502 (2005) [arXiv:astro- 
ph/0408426]. 

K. Jedamzik, Phys. Rev. D 74, 103509 (2006) [arXiv:hep-ph/0604251]. 

K. Jedamzik, K. Y. Choi, L. Roszkowski and R. Ruiz de Austri, JCAP 0607, 007 

(2006) [arXiv:hep-ph/05 12044]. 

M. Kusakabe, T. Kajino and G. J. Mathews, Phys. Rev. D 74, 023526 (2006) 
[arXiv:astro-ph/0605255] . 

M. Pospelov, Phys. Rev. Lett. 98, 231301 (2007) [arXiv:hep-ph/0605215]. 

S. Dimopoulos, D. Eichler, R. Esmailzadeh and G. D. Starkman, Phys. Rev. D 41, 
2388 (1990). 

A. De Rujula, S. L. Glashow and U. Sarid, Nucl. Phys. B 333, 173 (1990). 

J. Rafelski, M. Sawicki, M. Gajda and D. Harley, Phys. Rev. A 44, 4345 (1991). 

K. Kohri and F. Takayama, Phys. Rev. D 76, 063507 (2007) [arXiv:hep-ph/0605243]. 

M. Kaplinghat and A. Rajaraman, Phys. Rev. D 74, 103004 (2006) [arXiv:astro- 
ph/0606209]. 

R. H. Cyburt, J. R. Elhs, B. D. Fields, K. A. Olive and V. C. Spanos, JCAP 0611, 
014 (2006) [arXiv:astro-ph/0608562]. 

F. D. Steffen, AIP Conf. Proc. 903, 595 (2007) [arXiv:hep-ph/0611027]. 

J. Pradler and F. D. Steffen, Phys. Lett. B 648, 224 (2007) [arXiv:hep-ph/0612291]. 



26 



[20] K. Hamaguchi, T. Hatsuda, M. Kamimura, Y. Kino and T. T. Yanagida, Phys. Lett. 
B 650, 268 (2007) [arXiv:hep-ph/0702274]. 



[21 

[22 

[23 
[24; 

[25 
[26 
[27 

[28 
[29 
[30 



[31 
[32 
[33 
[34; 

[35 

[36 

[37 
[3 



C. Bird, K. Koopmans and M. Pospelov, arXiv:hep-ph/0703096. 

M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 649, 436 (2007) [arXiv:hep- 
ph/0703122]. 

F. Takayama, Phys. Rev. D 77, 116003 (2008) [arXiv:0704.2785 [hep-ph]]. 

T. Jittoh, K. Kohri, M. Koike, J. Sato, T. Shimomura and M. Yamanaka, Phys. Rev. 
D 76, 125023 (2007) [arXiv:0704.2914 [hep-ph]]. 

K. Jedamzik, Phys. Rev. D 77, 063524 (2008) [arXiv:0707.2070 [astro-ph]]. 

J. Pradler and F. D. Steffen, Phys. Lett. B 666, 181 (2008) [arXiv:0710.2213 [hep-ph]]. 

J. Kersten and K. Schmidt-Hoberg, JCAP 0801, Oil (2008) [arXiv:0710.4528 [hep- 
ph]]. 

J. Pradler and F. D. Steffen, Eur. Phys. J. C 56, 287 (2008) [arXiv:0710.4548 [hep-ph]]. 

K. Jedamzik, JCAP 0803, 008 (2008) [arXiv:0710.5153 [hep-ph]]. 

M. Kusakabe, T. Kajino, R. N. Boyd, T. Yoshida and G. J. Mathews, Phys. Rev. D 
76, 121302 (2007) [arXiv:0711.3854 [astro-ph]]. 

M. Kusakabe, T. Kajino, R. N. Boyd, T. Yoshida and G. J. Mathews, arXiv:0711.3858 
[astro-ph] . 

M. Pospelov, arXiv:0712.0647 [hep-ph]. 

M. Kawasaki, K. Kohri, T. Moroi and A. Yotsuyanagi, arXiv:0804.3745 [hep-ph]. 

F. D. Steffen, Phys. Lett. B 669, 74 (2008) [arXiv:0806.3266 [hep-ph]]. 

J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B 588, 7 (2004) 
[arXiv:hep-ph/0312262] . 

D. G. Cerdeno, K. Y. Choi, K. Jedamzik, L. Roszkowski and R. Ruiz de Austri, JCAP 
0606, 005 (2006) [arXiv:hep-ph/0509275]. 

K. Jedamzik, K. Y. Choi, L. Roszkowski and R. Ruiz de Austri, JCAP 0607, 007 

(2006) [arXiv:hep-ph/05 12044]. 

T. Asaka, K. Hamaguchi and K. Suzuki, Phys. Lett. B 490, 136 (2000) [arXiv:hep- 
ph/0005136]. 

M. Bolz, A. Brandenburg and W. Buchmiiller, Nucl. Phys. B 606, 518 (2001) [Erratum- 
ibid. B 790, 336 (2008)] [arXiv:hep-ph/0012052]. 

27 



[39 
[40 

[41 
[42 
[43 

[44; 
[45 

[46 

[47; 

[48 
[49 

[50 

[51 
[52 
[53 

[54; 



M. Fujii, M. Ibe and T. Yanagida, Phys. Lett. B 579, 6 (2004) [arXiv:hep-ph/0310142]. 

J. L. Feng, S. Su and F. Takayama, Phys. Rev. D 70, 075019 (2004) [arXiv:hep- 
ph/0404231]. 

F. D. Steffen, JCAP 0609, 001 (2006) [arXiv:hep-ph/0605306]. 

J. Pradler and F. D. Steffen, Phys. Rev. D 75, 023509 (2007) [arXiv:hep-ph/0608344]. 

V. S. Rychkov and A. Strumia, Phys. Rev. D 75, 075011 (2007) [arXiv:hep- 
ph/0701104]. 

F. D. Steffen, arXiv:071 1.1240 [hep-ph]. 

S. Ambrosanio, B. Mele, S. Petrarca, G. Polesello and A. Rimoldi, JHEP 0101, 014 
(2001) [arXiv:hep-ph/0010081]. 

J. R. Elhs, A. R. Raklev and O. K. Oye, JHEP 0610, 061 (2006) [arXiv:hep- 
ph/0607261]. 

W. BuchmiiUer, K. Hamaguchi, M. Ratz and T. Yanagida, Phys. Lett. B 588, 90 
(2004) [arXiv:hep-ph/0402179]. 

H. U. Martyn, Eur. Phys. J. C 48, 15 (2006) [arXiv:hep-ph/0605257]. 

K. Hamaguchi, M. M. Nojiri and A. de Roeck, JHEP 0703, 046 (2007) [arXiv:hep- 
ph/0612060]. 

M. Asplund, D. L. Lambert, P. E. Nissen, F. Primas and V. V. Smith, Astrophys. J. 
644, 229 (2006) [arXiv:astro-ph/0510636]. 

R. Cayrel et al, arXiv:0708.3819 [astro-ph]. 

F. Primas, M. Asplund, P. E. Nissen and V. Hill, arXiv:astro-ph/0009482. 

F. Primas, P. Molaro, P. Bonifacio and V. Hill, arXiv:astro-ph/0008402. 

A. M. Boesgaard and M. C. Novicki, Astrophys. J. 633, L125 (2005) [arXiv:astro- 
ph/0509483]. 

A. M. Boesgaard and M. C. Novicki, Astrophys. J. 641, 1122 (2006) [arXiv:astro- 
ph/0512317]. 



[55] A. J. Korn et al, Nature 442, 657 (2006) [arXiv:astro-ph/0608201]. 

[56] H. Utsunomiya et al, Phys. Rev. C63, 018801 (2000). 

K. Sumiyoshi, H. Utsunomiya, S. Goko, and T. Kajino, Nucl. Phys. A709, 467 (2002). 

[57] J. D. Garcia, N. H. Kwong, and J. S. Cohen, Phys. Rev. A 35, 4068 (1987). 

28 



[58] H. Reeves, W. A. Fowler, and F. Hoyle, Nature 226, 727 (1970). 

M. Meneguzzi, J. Audouze, and H. Reeves, Astron. Astrophys. 15, 337 (1971). 

[59] E. Vangioni-Flam, J. Audouze, Y. Oberto, and M. Casse, Astrophys. J. 364, 568 
(1990). 

[60] B. D. Fields, K. A. Olive, E. Vangioni-Flam and M. Casse, Astrophys. J. 540, 930 
(2000) [arXiv:astro-ph/9911320]. 

[61] B. D. Fields, K. A. Olive and E. Vangioni-Flam, Astrophys. J. 623, 1083 (2005) 
[arXiv:astro-ph/041 1728] . 

[62] V. F. Mukhanov, Int. J. Theor. Phys. 43, 669 (2004) [arXiv:astro-ph/0303073]. 

[63] G. R. Caughlan and W. A. Fowler, Atom. Data Nucl. Data Tabl. 40, 283 (1988). 

[64] E. Cremmer, S. Ferrara, L. Girardello and A. Van Proeyen, Nucl. Phys. B 212, 413 
(1983). 

[65] J. Wess and J. Bagger, Princeton, USA: Univ. Pr. (1992) 259 p 

[66] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[67] W. Buchmiiller, L. Covi, K. Hamaguchi, A. Ibarra and T. Yanagida, JHEP 0703, 037 
(2007) [arXiv:hep-ph/0702184]. 

[68] C. F. Berger, L. Covi, S. Kraml and F. Palorini, JCAP 0810, 005 (2008) 
[arXiv:0807.0211 [hep-ph]]. 

[69] W. Buchmiiller, K. Hamaguchi, M. Ibe and T. T. Yanagida, Phys. Lett. B 643, 124 
(2006) [arXiv:hep-ph/0605164]. 

[70] L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. Lett. 82, 4180 (1999) [arXiv:hep- 
ph/9905212]. 

L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, JHEP 0105, 033 (2001) [arXiv:hep- 
ph/0101009]. 

L. Covi, L. Roszkowski, R. Ruiz de Austri and M. Small, JHEP 0406, 003 (2004) 
[arXiv:hep-ph/0402240] . 

A. Brandenburg and F. D. Steffen, JCAP 0408, 008 (2004) [arXiv:hep-ph/0405158]. 

[71] A. Brandenburg, L. Covi, K. Hamaguchi, L. Roszkowski and F. D. Steffen, Phys. Lett. 
B 617, 99 (2005) [arXiv:hep-ph/0501287]. 

[72] K. Y. Choi, L. Roszkowski and R. Ruiz de Austri, JHEP 0804, 016 (2008) 
[arXiv:0710.3349 [hep-ph]]. 



29 



[73] K. Hamaguchi, Y. Kuno, T. Nakaya and M. M. Nojiri, Phys. Rev. D 70, 115007 (2004) 
[arXiv:hep-ph/0409248] . 

[74] J. L. Feng and B. T. Smith, Phys. Rev. D 71, 015004 (2005) [Erratum-ibid. D 71, 
019904 (2005)] [arXiv:hep-ph/0409278]. 

[75] M. Kamimura, Y. Kino and E. Hiyama, arXiv:0809.4772 [nucl-th]. 



30 



